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ABSTRACT 

We present new number density and surface brightness profiles for the globular cluster M92 
(NGC 6341). These profiles are calculated from optical images collected with the CCD mosaic camera 
MegaCam at the Canada-France-Hawaii- Telescope and with the Advanced Camera for Surveys on 
the Hubble Space Telescope. The ground-based data were supplemented with the Sloan Digital Sky 
Survey photometric catalog. Special care was taken to discriminate candidate cluster stars from field 
stars and to subtract the background contamination from both profiles. By examining the contour 
levels of the number density, we found that the stellar distribution becomes clumpy at radial distances 
larger than ~13', and there is no preferred orientation of contours in space. We performed detailed fits 
of King and Wilson models to the observed profiles. The best-fit models underestimate the number 
density inside the core radius. Wilson models better represent the observations, in particular in the 
outermost cluster regions: the good global agreement of these models with the observations suggests 
that there is no need to introduce an extra-tidal halo to explain the radial distribution of stars at large 
radial distances. The best-fit models for the number density and the surface brightness profiles are 
different, even though they are based on the same observations. Additional tests support the evidence 
that this fact reflects the difference in the radial distribution of the stellar tracers that determine the 
observed profiles (main sequence stars for the number density, bright evolved stars for the surface 
brightness) . 

Subject headings: globular clusters: general — globular clusters: individual (M92) — stars: kinematics 
and dynamics — stars: Population II 



1. INTRODUCTION 

The Galactic globular clusters (GCs) are the oldest 
(~ 11-13 Gyr) Galactic systems and have complex inter- 
nal and external dynamics (Gnedin & Ostriker 1997). 
The internal dynamics is driven by two-body relaxation 
with a time scale that is typically shorter than their age 
(Meylan & Heggie 1997). Therefore, the density profile 
of the innermost regions is expected to be well described 
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by King models (1962,1966). The outermost regions are 
characterized by the interaction with external tidal forces 
(Spitzer 1958; Spitzer & Chavalier 1973; Aguilar et al. 
1988) and by the evaporation of low-mass stars (Spitzer 
& Harm 1958). These phenomena are expected to pro- 
duce deviations from the spherical King models. By mov- 
ing beyond the cluster truncation radius (rt) the escaping 
stars can form halos or extended tidal tails. 

The first empirical evidence of extra-tidal structures 
in GCs was found from photographic plates (Grillmair 
et al. 1995; Lehmann & Scholz 1997; Testa et al. 2000 
[TOO]; Leon et al. 2000). More recently, investigations 
laased on CCD photometry found evidence either of tidal 
tails (Odenkirchen et al. 2001,2003; Grillmair & Johnson 
2006; Chun et al. 2010; Jordi & Grebcl 2010 [JGIO]) or 
of surrounding halos (Lee et al. 2003 [L03] ; Olszewski et 
al. 2009; JGIO; Correnti et al. 2011) around more than 
30 GCs. 

An accurate determination of the truncation radius of 
M92 and of the shape of its external regions is still miss- 
ing, even though several investigations have been car- 
ried out on this topic (see Table 1). The oldest Sur- 
face Brightness (SB) profile (obtained from photographic 
plates and small-format CCD images) was provided by 
Trager et al. (1995, [T95]). The same was later analyzed 
by McLaughlin & van der Marel (2005, [MLvdM05]). By 
using photographic plates, TOO found evidence of extra- 
tidal stars at 30' from the cluster center, and provided 
a surface density map, which shows only marginal ev- 
idence for an elongation orthogonal to the direction of 
the Galactic center. More recently, L03 by using a mo- 
saic CCD camera, confirmed the occurrence of extra- 
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tidal stars and showed that the marginal elongation ap- 
pears only for the brightest stars. JGIO analyzed the 
Sloan Digital Sky Survey (SDSS) photometric catalog 
and found the same elongated contours detected by TOO, 
even though their data do not cover the entire M92 area. 

2. PHOTOMETRIC DATASETS 

We used both ground-based data collected with the 36 
CCD mosaic camera MegaCam at the Canada-France- 
Hawaii- Telescope (CFHT) and space data collected with 
the Advanced Camera for Surveys (ACS) on the Hub- 
ble Space Telescope (HST). The MegaCam images were 
collected in the g\r\i',z' bandJ*^. To increase the radial 
extent of the sky area covered by our dataset, the ground- 
based data (MegaCam images) were supplemented with 
multiband {g,r,i,z) data collected by the SDSS (Aihara 
et al. 2011). As a whole the ground-based data cover 
~ 4° X 4° around the cluster center, but they do not 
uniformly cover the sky area around M92. 

The ACS data were collected with three different point- 
ings: pointing a. High Resolution Channel (HRC), six 
F435W (t=340 s each) and 155 F555W0 (with expo- 
sure times ranging from 10 to 100 s) images located 
across the cluster center; pointing f3, Wide Field Channel 
(WFC), three F606W (t=0.5,5,90 s) and three F814V^ 
(t=0.5,6,100 s) images located across the cluster center; 
pointing 7, WFC, three F475W (t=3,20,40 s) and three 
F814'WS (t=l, 10,20 s) images located 2' from the clus- 
ter center in the South-East direction. The reduction 
and the photometry of the F555W and F814W images 
have been presented and discussed by Di Cecco et al. 
(2010, [DiClO]). Panels a) and b) of Fig. □ display the 
sky coverage of the space- and ground-based datasets. 

To estimate the completeness of the CFHT data we 
adopted the data of the pointing 7 (ACS- WFC). We 
found that the completeness for i <22 mag is 54% for 
75<R<150 arcsec, 82% for 150<R<200 arcsec, and com- 
plete for larger distances. We applied the above com- 
pleteness corrections to the star counts (DiClO). We also 
estimated the completeness of the SDSS data by using 
the CFHT data and we found that for radial distances 
larger than 700 arcsec and i < 22 mag they are complete. 
For radial distances smaller than 75" we adopted the data 
of pointing /3. These data display a gap of 2.5" between 
the two CCDs. To fill this gap we selected two regions of 
2.5" at the edges of the gap and we randomly extracted 
half of the stars in each of the two regions. Once the gap 
was filled, we assumed that the data of pointing /3 located 
between 20" and 75" are complete (DiClO). The com- 
pleteness of pointing /3, for radial distance -R- smaller 
than 11", was estimated using pointing a. In this case 
we found that in the magnitude range between the Main 
Sequence Turn Off (MSTO) and i^5551^-21.5 mag the 
pointing (3 dataset was complete at the 70% level. The 
star counts in this cluster region were corrected account- 
ing for the above completeness correction. 

In summary, we are dealing with three different 

1)- ACS-WFC, pointing /3, -R < 1.25'. For radial 
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distances smaller than 0.18' the completeness was esti- 
mated using ACS~HRC, pointing a. For radial distances 
0.18' <R< 1.25' we only filled the gaps. 

2y CFHT, 1.25' < R < 30'. For radial distances 
1.25' < R < 3.33' the completeness was estimated using 
ACS-WFC, pointing 7. The comparison of CFHT with 
pointing 7 indicates that the former dataset for i < 22 
mag is complete at larger distances. 

SDSS, 30' < i? < 2°. The completeness was esti- 
mated using CFHT data; they are complete for i < 22 
mag and radial distances larger than 700 arcsec. How- 
ever, they do not uniformly cover the sky region around 
the cluster (see panel b) in Fig. 1) and the star counts 
were accordingly corrected. 

Panels a), b) and c) of Fig.[2]show the color-magnitude 
diagrams (CMDs) of pointing a, /3 and 7. Stars plot- 
ted in these panels were selected according to the sharp- 
ness (|sh|<l) DAOPHOT index (see, e.g.. Stetson 1987, 
1994). These CMDs display well defined sequences in 
the evolved phases as well as along the Main Sequence 
(MS). The Red Giant Branch (RGB) stars brighter than 
red Horizontal Branch (HB) stars {F555W^1A mag) are 
saturated in pointing a. 

To provide homogeneous star counts across the entire 
GC, the i- and the r-band from the SDSS, as well as the 
F8I4W- and the F606W-hand of pointing (3 were trans- 
formed into the i'- and the r-band of the MegaCam pho- 
tometric system. The accuracy of the transformations 
is better than 0.02 mag (DiClO). In the following, we 
use the g,r,i,z bands (without prime) to refer to the 
CFHT bands. The reason for the above transformations 
is threefold: a) i- and r-band are common to the different 
datasets; b) the data in these bands have good photomet- 
ric accuracy (CTi_i=0.06 mag, at least three magnitudes 
fainter than the MSTO); c) the i-band is minimally af- 
fected by saturation problems and it was adopted to com- 
pute the density profiles. 

3. RADIAL DENSITY PROFILE 

The field of view of the CFHT dataset (1° x 1°) and 
of the SDSS dataset fully enclose the estimated radial 
extent of M92 (see Table 1). In order to constrain the 
cluster edges, candidate cluster and field stars must be 
distinguished. The method used to identify the two dif- 
ferent groups of star is described below. 

We selected the photometric catalog by using the in- 
trinsic photometric error (a^-i <0.10 mag), the separa- 
tior(3 (sep>2.5), and the distance from the cluster center 
(10" < R <180"). We computed a fiducial line (ridge- 
line) in the z,r-i CMD by means of a three-dimensional 
Hess diagram (Ferraro et al. 2013, in preparation). Note 
that the above selection criteria were only applied to es- 
timate the ridge line. 

Panel a) of Fig. [3] shows the i,r-i CMD for the en- 
tire sample of stars together with the above ridgeline. 
Stars plotted in this figure were selected according to 
the radial distance and to a very mild photometric er- 
ror cut (cTi— i <0.25 mag). Note that in this analysis 
we only adopted space and CFHT data; the SDSS data 
will be discussed in the following. The solid line shows 

1^ The separation index quantifies the degree of crowding, i.e. 
the amount of spurious light, due to neighboring stars, affecting 
the magnitude of individual stars (Stetson et al. 2003). 
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the computed ridgeline, whereas the dashed ones display 
the acceptance region, that is, the region of the CMD 
where we assume that the candidate MS and RGB clus- 
ter stars lie. The acceptance region is centered on the 
ridgeline and the width in color goes from 0.01 mag close 
to the tip of the RGB up to 0.30 mag for magnitudes 
fainter than the MSTO. We cut the acceptance region at 
magnitudes brighter than «=21.7 mag, because at fainter 
magnitudes the photometric error in the color increases 
and the ridgeline is less well determined. 

Panel b) of Fig. [3] shows the accepted M92 stars, that 
is, the candidate MS, RGB, and HB stars. The MS and 
RGB accepted cluster stars include some candidate field 
stars with colors and magnitudes similar to those of M92. 
The HB stars, instead, can be easily selected since they 
are bluer than the field stars. Panel c) of Fig. |3] shows 
the GMD of candidate field stars (rejected) with their 
typical peaks in color around r-i=0-0.2 and r-i=1.2-1.4 
mag. The stars located at (i,r-i)'^(17.5,0.0) are M92 
Blue Stragglers. They were not included among accepted 
stars, because the fainter ones partially overlap with field 
and MSTO stars in the CMD, while the brighter ones are 
a minimal fraction of cluster stars. 

To further constrain the radial extent of candidate clus- 
ter stars we investigated the ratio between the number 
of accepted stars and the total number of stars in ra- 
dial bins located in the outermost cluster regions. Data 
plotted in Fig. |4] show the i,r-i CMD for stars in ra- 
dial bins located between 13' and 2° from the cluster 
center. The first three radial bins -panel a), b) and 
c- are entirely included inside the CFHT dataset, while 
the last two -panel d) and e)- are entirely located in- 
side the SDSS dataset. The two solid lines display the 
acceptance region we defined in Fig. [3l The photomet- 
ric precision in the outermost cluster regions is clearly 
supported by the thin distribution of MS stars. Data 
plotted in the three innermost radial bins indicate that 
MS stars are crucial to trace the radial extent of can- 
didate cluster stars. Moreover, the ratio between the 
number of cluster stars and the total number of stars 
is steadily decreasing when moving toward the outer- 
most cluster regions. It decreases from almost 50% 
{Na/Nt=0A9±0.02) for R~ 14.5' to slightly less than 
one third {Na/Nt=0.30±0.01) for R- 25'. This radial 
distance appears to be a preliminary plausible lower limit 
for the truncation radius, and indeed the same ratio in 
the two outermost radial bins attains smaller constant 
values. Note that the main vertical sequence partially 
overlapping with the acceptance region is almost entirely 
made up of field stars. 

To further constrain the plausibility of the above work- 
ing hypothesis concerning the radial extent of M92, we 
decided to investigate the radial distribution of extra- 
galactic sources. We adopted the entire set of r-band im- 
ages collected with MegaCam at CFHT and performed a 
new independent photometry by using Sextractor (Bertin 
& Arnouts 1996). The non-point like sources were se- 
lected -following Evans et al. (2010)- as the objects 
with a local point spread function 90% enclosed counts 
fraction larger than 1.4". The adopted value was fixed 
by eye inspection of the mean r-band image. The Fig. [S] 
shows the radial distribution of the non-point like sources 
(blue dots) over the mean r-band image with apparent 
magnitudes between r'--^19 and r~22 mag. To help the 



eye to identify the sky area covered by M92, the red and 
the purple circles display the Wilson truncation radius 
based on both the number density and the surface bright- 
ness profile (see §4). The smaller orange circles show the 
candidate galaxy clusters identified from the SDSS DR6 
(Wen et al. 2009). A glance at data plotted in this fig- 
ure shows that the candidate galaxy clusters are located 
either close to the truncation radius or beyond it. The 
innermost candidate galaxy cluster appears a bit suspi- 
cious, since it is located in a cluster region with a high 
stellar density. 

We selected the objects located inside the sky region 
covered by the five candidate galaxy clusters beyond the 
truncation radius and plotted them in the CMD and we 
found that a significant fraction (~65%) of them are lo- 
cated outside the acceptance region (blue triangles in 
panel c) of Fig. |4]). This fraction agrees, within the er- 
rors, quite well with the ratio between accepted and total 
number of stars (see the discussion above in this section). 

To remove spurious stars that were erroneously ac- 
cepted, we used the method described by Walker et al. 
(2011) for the GC IC 4499. This method works partic- 
ularly well for our datasets, due to the large sky area 
they cover. To estimate the density of the rejected stars, 
the sky area covered by CFHT (i? <30') and by SDSS 
(0.5°< R <2°) data was divided into concentric annuli. 
The star counts based on SDSS data were corrected to 
account for the non homogeneous coverage of this photo- 
metric catalog. Panel a) of Fig. [6] shows the logarithmic 
surface density of these objects (number of rejected stars, 
Nji, per arcmin^) as a function of the inverse of the ra- 
dial distance. We performed a linear fit to the individual 
points and by extrapolating to infinite radial distance we 
found that the asymptotic value is /ii =0.42±0.10 [loga- 
rithmic number of stars per arcmin^]. We also estimated 
the asymptotic value as the mean of the five outermost 
values, finding /i2=0.36±0.10. We adopted the mean of 
the above estimates: /^=0.39±0.14. 

In panel b) of Fig. [6] we plotted the logarithm of the 
ratio between the number of accepted stars and the num- 
ber of rejected stars {Na/Nh). This ratio was estimated 
using the entire dataset. To estimate the mean asymp- 
totic value, we adopted the average of the outermost 
three radial bins, obtaining x=~0.41±0.03. Eventually, 
by multiplying the number of rejected stars per arcmin^ 
by Na/Nji, we found the number of candidate field stars 
that were erroneously classified as candidate M92 stars. 
Notably, we found that the asymptotic number of spu- 
riously accepted stars is 10'' x 10-^ ~ 0.95 star/arcmin^. 
By subtracting this value from the number of the ac- 
cepted stars per unit area, we obtained the final Count 
Catalog of the candidate M92 stars. This catalog was 
used to compute the Number Density (ND) radial pro- 
file. We divided the cluster into concentric annuli and 
we counted the number of stars per arcmin^ that fall in- 
side each region, obtaining a radial profile ranging from 
R ^ 1.5" out to R ^2°. The error on each of the points 
is calculated as the square root of the number of stars, 
divided by the area of the annulus. 

Following Walker et al. (2011), we measured the spu- 
rious stellar flux which affects the stellar luminosity of 
accepted stars. The logarithmic flux density of the re- 
jected stars per arcmin^ (FIuxr) and the ratio between 
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the flux of accepted and rejected stars {Flux a/ Flux ji) 
are plotted in panels c) and d) of Fig. [6] The logarithm of 
the surface flux density of the rejected stars approaches 
ei=-5.72±0.03 [logarithmic star flux per arcmin^] when 
extrapolated to infinite radial distance, whereas it is 
e2=-5.76±0.03 when we use the mean value of the last 
three points. The mean of these values is e=-5.74±0.04, 
whereas the logarithm of the ratio between accepted 
and rejected stars is ^=-0.82±0.07. In this case, the 
flux of spurious accepted stars is 10^ x 10^ ~ 0.3 10~^ 
Flux/arcmin^. By subtracting this last value from the 
accepted stellar flux density, we obtained an independent 
flnal Flux Catalog for the candidate M92 stars. This 
catalog was used to calculate the SB radial profile. As 
in the previous case we divided the cluster into annular 
regions, and the SB for each annulus was computed by 
adding the flux contribution of the stars located inside 
the annulus and by dividing for its area. Data plotted in 
panels b) and d) of Fig. [5] show that the ratio between 
the number of accepted and rejected stars is more robust 
than the ratio between the flux of accepted and rejected 
stars, since the intrinsic dispersion of the former one is 
at least a factor of two smaller than the latter one. The 
difference is caused by the fact that the number ratio is 
rooted in the radial distribution of MS stars, while the 
flux ratio traces the radial distribution of bright evolved 
stars. To further constrain the role of bright stars in de- 
termining the radial slope of the SB, we calculated two 
more SB profiles by considering only stars fainter than 
a limiting magnitude of (SB-15) and i—17 (SB-17) 
mag (see Sect. S]). 

The error on the surface brightness profile was esti- 
mated by propagating the error on individual measure- 
ments of star magnitudes. The intrinsic photometric er- 
ror, in the magnitude range adopted to estimate the sur- 
face brightness profile, is typically of the order of a few 
hundredths of a magnitude (see § 2), because ground- 
based images were collected in good seeing conditions 
(DiClO) and we typically have more than ten images 
per band. The same applies for ACS images adopted in 
the central regions. We also calculated the error of the 
absolute photometric zero-points. Following DiClO, we 
calibrated the CFHT photometric catalog by using the 
local standards by Clem et al. (2007). The ACS and the 
SDSS photometric catalogs were also transformed into 
the same photometric systems by using the new local 
standards. We ended up with a mean calibration error 
in the i-hand of 0.02±0.04 mag (Di Cecco 2009). This 
error was eventually summed in quadrature with the in- 
trinsic photometric error. 

We evaluated the symmetry of the ND as a function of 
the radial distance by using the Count Catalog. To avoid 
possible systematic uncertainties in the radial distribu- 
tion, the symmetry of the ND was estimated on the basis 
of ACS and CFHT data. The SDSS data were neglected 
since they do not uniformly cover the area of the sky 
around the cluster center. The conclusions concerning 
the departure from circular symmetry of the radial dis- 
tribution for distances between 6' and 30' are not affected 
by the inclusion of the SDSS dataset. We computed the 
contour levels of the candidate M92 and field stars (black 
and red lines in panel a) of Fig. [7]). The contour levels 
become less circular symmetric when moving toward the 
outermost regions, and for radial distances between 6' 



and 10' the stellar density decreases by almost one order 
of magnitude (see top panels of Fig. [7]). The contour lev- 
els become asymmetric at a distance of ~13' (see panel 
a) of Fig. [T]) , which is almost equivalent to the truncation 
radius of M92 (see column 5 in Table 1) available in the 
literature. Outside this region the distribution of candi- 
date cluster stars becomes clumpy. Panels b) and c) of 
Fig.[7]show the projected linear density (marginal) along 
the axes of the candidate M92 and field stars (black and 
red lines). By inspecting these panels, we notice that 
the detected sharp decrease in density is associated with 
the radial distance at which the surface density of the 
candidate M92 stars becomes smaller than that of the 
candidate field stars. 

To trace in detail the departure of the contour levels 
from circular symmetry, we performed a fit of each con- 
tour with a circle. The center of the circles is identical to 
the center of the contours -i.e. the cluster center- and 
the radius of the circles is the fitting parameter. Then, 
we computed the residuals in arcsec between individual 
contours and best fitting circles. The residuals were esti- 
mated from the very center of the cluster out to a radial 
distance of i? ~ 800". 

The residuals plotted in panel d) of Fig.[7]show that the 
innermost contour levels display symmetric radial distri- 
butions, and indeed the residuals attain vanishing values 
out to R ^3'. At larger radial distances, the asymme- 
try increases out to R -9'-10' {R -500-600"), where 
the residuals show a shoulder clearly connected with the 
density drop detected in the contour plot. At even larger 
distances the contours become more asymmetric (resid- 
uals ~ 26" for distances of i? ~ 13'); the fit in the outer- 
most regions fails to converge because of the large asym- 
metries. The increase in asymmetry that we found in the 
outer regions could be the consequence of the ffuctuations 
associated with the decrease in density. To validate this 
working hypothesis we performed a series of simulations 
by using the observed density profile to compute syn- 
thetic GCs; the radial distribution of the synthetic GCs 
was required to be symmetric. We applied to these GCs 
the same procedure to evaluate the contour levels and 
the same fit with circles of variable radius. The vertical 
hatched area plotted in panel d) of Fig. [7] marks the resid- 
uals calculated for the synthetic GCs. The comparison 
between this area and the plotted points indicates that 
the asymmetry in the real cluster is at least 3cr larger 
than in the synthetic clusters. To characterize further 
the nature of the asymmetries in the contour levels, we 
show two green arrows, in the top panel of Fig. [71 indi- 
cating the direction of the Galactic center (long arrow) 
and of the M92 proper motion (short arrow) according 
to Dinescu et al. (1999). We found no clear correlation 
between these directions and the clumpy distribution of 
candidate M92 stars at large radial distances. 

Evidence of a clumpy stellar distribution in the out- 
skirts of M92 was also present in the stellar density maps 
provided by Testa et al. (2000, see their Fig. 6), by Lee 
et al. (2003, see their Fig. 12) and by Jordi & Grebel 
(2010, see their Fig. 17). The above results support the 
findings by Lee et al. (2003), concerning the marginal 
evidence of an elongation of outermost clumpy stars in 
the direction orthogonal to the direction of the Galactic 
center. 
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4. DYNAMICAL MODELS AND FITS 

We carried out fits of dynamical models to the observed 
radial profiles. We considered the King (1966) and the 
Wilson (1975) spherical and isotropic dynamical models, 
defined by the following distribution functions: 

Jk^A (e-"^ - e--^^") E<Eo, (1) 

/w = ^{e-'^^-e-"^"[l-«(^-£^o)]} E<Eo. (2) 

The quantity E is the specific star energy E — + 
$(r), where $(r) is the mean-field gravitational poten- 
tial, to be determined from the Poisson equation. Both 
distribution functions vanish for energies larger than the 
threshold energy Eq, corresponding to stars to be consi- 
dered as unbound. The energy truncation can be trans- 
lated into a truncation radius, rt, which indicates the 
boundary of the system. For each family of models, the 
constants A, Eq, and a in Eqs. (P) and (0) define two 
dimensional scales (a typical radius and a typical mass 
or velocity) and one dimensionless parameter, the central 
depth of the potential well (related to the concentration 
parameter). We recall that fitting by a one-component 
dynamical model assumes that the underline stellar pop- 
ulations are distributed homogeneously. To identify the 
best-fit model we adopted the procedure described by 
Zocchi et al. (2012). The results are shown in Figs. [8] 
and m the values of the relevant parameters of the best- 
fit models are listed in the upper part of Table [H 

The fit to the ND profile is shown in the top panel of 
Fig. m Both King and Wilson models underestimate the 
central ND profile, failing with respect to the four inner- 
most points. A quantitative interpretation of this dis- 
crepancy remains unavailable, but the problem is likely 
to be related to the failure of the assumptions at the basis 
of a one-component dynamical model in the central re- 
gions. Concerning the King models, the good agreement 
with observations that can be found in the middle part 
of the profile breaks down around 700", where the pro- 
file approaches the background level (0.95 star/arcsec^). 
In the case of the Wilson models, instead, only the two 
outermost points are discrepant, and the model fits the 
data out to a distance greater than 1000". Note that 
the two outermost points are likely to be affected by er- 
rors in the subtraction of background stars. The bottom 
panel of Fig. |8] shows the fits to the SB profile. When 
compared to the King best-fit model, the Wilson best- 
fit model provides a more adequate overall description, 
not only in relation to the outermost points, as in the 
ND profile, but also in the central part of the profile. 
The satisfactory performance of the Wilson models indi- 
cates that the observations can be explained by means 
of a less abrupt truncation radius, with no need to intro- 
duce extra-tidal halos. In this case, the two outermost 
points were not taken into account to calculate the best- 
fit parameters, since they are expected to be even more 
affected by errors in the subtraction of background stars, 
compared to the corresponding points in the ND profile. 
In the figure, the background level (25.30 mag/arcsec^) 
is indicated as a horizontal dashed line. 

Surprisingly, even if the ND and the SB profiles come 
from the same set of observations, the best-fit parameters 
determined by the fits are significantly different. We ar- 
gue that this behavior is due to the fact that each profile 



represents a different aspect of the density distribution 
of the cluster. On the one hand, the ND profile, derived 
by considering the radial distribution of both luminous 
and faint stars, is dominated by the MS stars, which 
greatly outnumber evolved (RGB, HB) stars (Castellani 
et al. 2007). On the other hand, the SB profile is heav- 
ily affected by the presence of the brighter RGB stars. 
The difference in the best-fit models reflects the intrinsic 
difference in the radial distribution of the stellar tracer 
that determines each profile. This behavior should be 
interpreted as a signature of mass segregatioiO- Indeed, 
the evolved slightly more massive star^ij appear to be 
more centrally concentrated compared to stars with lower 
masses (see the values of the concentration parameter c 
in Table H]). 

This interpretation is confirmed by an additional test. 
We carried out the same fitting procedure on the SB- 15 
and the SB-17 profiles (as defined in Sect. 3). The result- 
ing parameters are listed in the third and fourth rows of 
Table [T] The fits are shown in the top and in the bot- 
tom panel of Fig. |9l respectively. Inspecting the values 
of the best-fit parameters, it appears that by eliminat- 
ing the brightest stars, the profiles tend to approach the 
ND profile. Indeed, the values of the concentration pa- 
rameter and of the radial scale -r^- follow a monotonic 
trend from the SB profile to the SB-15, to the SB-17 and 
finally to the ND profile. Interestingly enough, for the 
SB-15 and SB-17 profiles, the King models reproduce the 
data better than the Wilson models, in contrast to what 
we found for the previously described ND and SB pro- 
files. A reason for this can be found by comparing the 
radial extent of the different profiles. The three surface 
brightness profiles are shaped by the fact that the cut of 
candidate cluster stars brighter than a limiting magni- 
tude causes a decrease in the radial extent and a more 
abrupt truncation of the profile. Indeed, the outermost 
radial point for the SB-15 profile is located at R ^ 700", 
and for the SB-17 profile is located at i? - 500". The 
larger radial extent of the SB profile is the consequence 
of a few bright giants that keep the SB profile well above 
the background level. 

For comparison, the middle part of Table [1] lists the 
most recent best-fit parameters published for M92. Note 
that in only a few cases the Wilson parameters are avail- 
able. The acronyms in the first column identify the pa- 
pers in which the results have been published. 

The lower part of Table [1] lists the results of the fits 
to the other twcEB available surface brightness profiles in 
the literature, the one by T95, and the one by LOS. We 
decided to add to both datasets the surface brightness 
data by Noyola et al. 2006 (N06), which are relative to 
the innermost region of the cluster. The results of these 
fits are shown in Fig. [101 For the T95+N06 profile, the 
models provide an equally good fit to the data; we note 
that in this case the outermost point in the profile is 
located at a radial distance R< 700". For the L03-HN06 
profile, the good agreement with observations that is seen 
in the inner and middle parts of the profile breaks down 

The occurrence of mass segregation in M92 was also suggested 
by Andreuzzi et al. (2000). 

Note that hot HB stars are less massive than MSTO stars, 
but they are a minor fraction of evolved cluster stars. 

Another surface brightness profile was presented by JGIO, but 
the dataset is not publically available. 
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around 600" for the King model, and at more than 1000" 
for the Wilson model. We recall that a similar behavior 
was also found for the fit to the ND profile previously 
described. In the last profile the outermost points are 
not reproduced by any of the considered models. 

The values of the parameters listed in Table [T] are not 
consistent with each other, within the errors. In conclu- 
sion, we believe that a proper comparison of the values 
of the best-fit parameters found by fitting models to dif- 
ferent (ND, SB) profiles requires that we take into ac- 
count the role of the different stellar tracers in determin- 
ing their shape, which makes standard one-component 
dynamical models questionable. 

5. DISCUSSION AND FINAL REMARKS 

We studied the radial distribution of stars of the glob- 
ular cluster M92 by using ground-based (MegaCam at 
CFHT, SDSS) and space (ACS on HST) data. 

The contour levels, based on star count data, are sym- 
metric in the innermost regions, and exhibit an increas- 
ing asymmetry for radial distances between 3' and 9'- 
10'. For distances larger than ~13' the stellar distribu- 
tion becomes clumpy. The contour levels do not exhibit 
a preferred orientation in space. We calculated two in- 
dependent radial profiles, to describe the distribution of 
stars in the cluster, the number density (ND) and the 
surface brightness (SB) profile. To calculate these pro- 
files, we subtracted the background contamination with 
two independent methods. We performed fits of spherical 
King and Wilson models to the above profiles. Wilson 
models appear to reproduce, better than King models, 
the behavior of the outermost regions of the cluster with 
no need of extra-tidal halos. Interestingly, for the ND 
profile, both models significantly underestimate the ob- 
servations in the innermost regions. 

We also found that the best fit to the ND and to the SB 
profile are provided by two different models for the two 
families, even though the profiles are derived from the 
same datasets. We argue that this difference is caused by 
a difference in the radial distribution of the stellar trac- 
ers that characterize the two observed profiles. The ND 
profile traces the radial distribution of MS stars, whereas 
the SB profile that of bright evolved (RGB, HB) stars. 
This conclusion is supported also by the results of a test 
that has been carried out on two additional profiles, cal- 
culated by considering only stars fainter than a given 



magnitude. 

Hopefully, a thorough discussion of the behavior in this 
and other clusters should determine which of the various 
profiles considered is best suited for a study in terms 
of one-component models. To our knowledge, this is the 
first investigation in which independent estimates of both 
ND and SB profiles are provided, starting from the same 
set of data, and compared. 

In this context a key role can be played by the new 
generation of wide field imagers that are available at the 
4-8m class telescopes (Dark Energy Camera and Survey 
at the CTIO 4m Blanco telescope, Mohr et al. 2012; 
Hyper SuprimeCam at SUBARlH). In a single pointing 
they can cover the entire extent of a large number of 
GCs and with modest exposure times will allow us to 
perform homogeneous and accurate photometry several 
magnitudes fainter than the MSTO. 
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TABLE 1 
Best-fit parameters. 



Profile 




King 


models 






Spherical Wilson models 








c 


rs 


rt 




c 


rs 


n 


ND 


6.91±0.02 


1.50±0.01 


34.25±0.38 


18.11±0.44 


5.84±0.02 


1.73±0.01 


48.03±0.42 


42.63±1.19 




o.40±0.Ui 


i.y5±u.ou 


1 c oo_l_r> no 
iO.zzitU.Uz 


no Qn_l_n i o 


O.DO±O.U1 


2.14±0.00 


21.08±0.05 


48.8U±U.80 


SB-15 


7.20±0.01 


1.59±0.00 


19.43±0.05 


12.56±0.10 


6.52±0.01 


2.06±0.00 


24.33±0.06 


46.67±0.42 


SB-17 


6.95±0.01 


1.51±0.00 


23.17±0.06 


12.63±0.12 


6.29±0.01 


1.94±0.00 


29.04±0.06 


41.97±0.53 


T95 


7.92 


1.81 


23.67 


15.20 










TOO 








12.33 










LOS 


8 


1.83 


12.42 


14.00 










MLvdM05 


7.5 


1.68 


16.15 


12.88 


5.9 


1.75 


26.51 


24.85 


JGIO 


6.93 


1.51 


23.37 


12.55 










T95+N06 


7.54 


1.69 


14.72 


12.07 


6.34 


1.96 


19.33 


29.58 


L03-I-N06 


7.84 


1.78 


13.46 


13.67 


6.61 


2.12 


18.94 


41.72 



Note. — For each model, we list the dimensionless parameter ^f, the concentration c, the scale radius rs (arcsec), and the 
truncation radius rt (arcmin). The different profiles are identified by the label in the first column. The cases indicated as 
T95+N06 and L03+N06 refer to the fits we performed on composite profiles, obtained by combining the profiles from T95 and 
LOS with the profile from Noyola et al. (2006, [N06]), which covers the innermost regions of the cluster. 



9 



43.4 



43.5 



^' 43.2 



§ 43.1 



^5 



43.0 



42.9 



1 1 1 n n I 1 1 i-i 1 1 1 1 I 1 1 I I 1 1 1 1 ! I 1 1-| 1 1 I.1 1 I I I.I 1 1 I 1 1 I 1 1 I 1 1 1 I i-i I I.I .1 1 



HRC ft 



WFC j? 
WFC > 




I -I I n I , I r, 



■ II- I'l 



259,5 259,3 259,1 
a (J2000) [deg] 




III I I 



rf^m.-jyi i., ;-|,uj.'^i I I ,1. i.fluliS: M if-L i-u Hi ,T,1, 



262 261 260 259 258 257 
a {J2000) [deg] 



Fig. 1. — Panel a) - area of the sky across the globular cluster M92 covered by the different sets of space images collected with the 
Advanced Camera for Surveys (ACS) on board the HST. The field of view of the sky plot is 1° X 1°. The red and the blue squares show 
the images collected with the Wide Field Channel (WFC, pointings 13, 7), while the green square those collected with the High Resolution 
Channel (HRC, a). The black dots display the photometric catalog based on ground-based images collected with CFHT. The orientation 
is shown in the bottom right corner. Panel b) — same as Panel a), but for datasets collected with ground-based telescopes, namely CFHT 
(black dots) and SDSS (brown dots). Note that the latter dataset does not uniformly cover the area of the sky around M92. The field of 
view of the sky plot is 4° X 4° . See text for more details. 
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Fig. 2. — CMDs of M92 based on ACS at HST, pointings a, /3, and 7 in panels a), b), and c), respectively. Stars plotted in these CMDs 
were selected according to different selection criteria. Data of pointing 7 are the same as in DiClO. 




Fig. 3.— Panel a) i,r-i CMD based on both ground-based (MegaCam at CFHT) and space (ACS-WFC at HST) data. Panels b) and c): 
CMDs for candidate M92 {accepted) and candidate field (rejected) stars, respectively. See text for more details. 
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N^/Nt^=0.49 ±0.02 Na/Nt^=0.3 6±0.02 

i3'<R<ie' 1 e'<R<20' 

o) // b) /; 



O-O 0.5 




Fig. 4. — Panel a) - i,r-i CMD based on MegaCam/CFHT data. The red solid lines display the acceptance region defined in Fig. [3] 
The radial extent of the bin and the ratio between candidate cluster stars and total number of stars are also labeled. Panels b) and c) - 
Same as for panel a), but for stars located at larger radial distances. The blue triangles plotted in panel c) are the objects located inside 
the candidate galaxy clusters (Wen et al. 2009). See text for more details. Panels d) and e) - same as for panel a), but for stars entirely 
located inside the SDSS dataset. 
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Fig. 5. — M92 r-band image of the CFHT dataset. The blue dots display the non-point like sources identified using Sextractor. The 
large purple and red circles display the Wilson truncation radius according to the SB and the ND profile. The small orange circles show 
the candidate galaxy clusters identified by Wen et al. (2009). The orientation is: North up and East left, the field of view is 1° X 1°. 
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Fig. 6. — Panels a) and b): logarithm of the number density of the rejected stars and of the ratio between the number of accepted and 
rejected -Nj^/Nn— stars as a function of the inverse of the radial distance. Panels c) and d): logarithm of the flux density of the rejected 
stars, and of the ratio between the flux of accepted and rejected — Flux j^/ Flux ji- stars. The vertical arrows display the edge of outermost 
annulus, the solid lines show the linear fits over the entire samples, the dashed lines indicate the mean of the outermost annuli (red points). 
See text for more details. 
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Fig. 7. — Panel a): contour levels for candidate M92 (black) and field (red) stars projected onto the sky. The long green arrow marks 
the direction of the Galactic Center, the short one the M92 proper motion. Panels b) and c): projected logarithmic distribution (marginal) 
along the horizontal and the vertical axis for candidate M92 (black) and field (red) stars. Panel d): the black filled circles show the residuals 
of the fits to the contour levels plotted in panel a) with circles of variable radius as a function of the radial distance. The vertically hatched 
area shows the results of simulations (see text for more details). 
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Fig. 8. — Fits by King and Wilson spherical models to the ND (top) and SB (bottom) profiles. Solid lines correspond to the King-model 
fits, dotted lines to Wilson-model fits; the horizontal dashed line shows the background level; errors are shown as vertical error bars. Red 
squares indicate data obtained from pointing 7, grey triangles those from the CFHT dataset, and green circles those from the SDSS. Empty 
symbols mark regions for which a completeness correction was applied. 
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Fig. 9. — Fits by King and Wilson spherical models to the SB-15 (top) and the SB-17 (bottom) profiles. Solid lines correspond to the 
King-model fits, dotted lines to Wilson-model fits; the horizontal dashed line shows the background level. Data-points are indicated as in 
Fig.[8l and errors are shown as vertical error bars. 
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Fig. 10. — Fits by King and Wilson spherical models to the composite T95+N06 (top) and to the L03+N06 (bottom) profiles. Solid lines 
correspond to the King-model fits, dotted lines to Wilson-model fits; errors are shown as vertical error bars. Blue triangles indicate data 
from N06, red crosses those from other sources. 



